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Abstract. A general formalism to evaluate time-reversal odd transverse momentum 
dependent parton distributions (T-odd TMDS) is reviewed and applied to the calculation of 
the leading-twist quantities, i.e., the Sivers and the Boer-Mulders functions. Two different 
models of the proton structure, namely a non relativistic constituent quark model and the MIT 
bag model, have been used. The results obtained in both frameworks fulfill the Burkardt sum 
rule to a large extent. The calculation of nuclear effects in the extraction of neutron single 
spin asymmetries in semi-inclusive deep inelastic scattering off transversely polarized He is 
also illustrated. In the kinematics of JLab, it is found that the nuclear effects described by an 
Impulse Approximation approach are under theoretical control. 



1. Introduction 

The partonic structure of transversely polarized hadrons is still rather poorly known pQ. As a 
matter of facts, despite its leading-twist character, the transversity distribution is chiral-odd in 
nature and not accessible therefore in inclusive esperiments. Due to going-on and forthcoming 
important measurements in various Laboratories, this issue is nevertheless one of the most widely 
studied by the hadronic Physics Community, in particular by the Italian one, which is providing 
crucial contributions to a better understanding of this fascinating subject (see, for example, 
Refs. j2| and references therein). 

Semi-inclusive deep inelastic scattering (SIDIS) is one of the proposed processes to access 
the parton structure of transversely polarized hadrons. The theoretical description of semi- 
inclusive processes implies a more complicated formalism, accounting for the transverse motion 
of the quarks in the target O HI [5]. In particular, the non-perturbative effects of the intrinsic 
transverse momentum kx of the quarks inside the nucleon may induce significant hadron 
azimuthal asymmetries [6j [7] . 

The Sivers and the Boer-Mulders functions were defined in this scenario [U |9]. Transverse 
Momentum Dependent pdfs (TMDs) are the set of functions that depend on the intrinsic 
transverse momentum of the quark, in addition to the dependences, typical of the PDs, on the 
Bjorken variable and on the momentum transfer Q 2 . Their number is fixed counting the scalar 
quantities allowed by hermiticity, parity and time-reversal invariance. However, the existence of 
final state interactions (FSI) allows for time-reversal odd functions |10j . In effect, by relaxing 
this constraint, one defines two additional functions, namely, the Sivers and the Boer-Mulders 
(BM) functions. In SIDIS, the Sivers function is involved in the description of the single spin 
asymmetry, measured when an unpolarized beam scatters transeversely polarized targets. The 



single spin asymmetries are obtained constructing the difference of semi-inclusive cross-sections 
with different transverse polarization of the target with respect to the momentum transfer. The 
BM function appears in the azimuthal asymmetry in unpolarized SIDIS. The latter object refers 
to the detection of the produced hadron at different angles with respect to the plane containing 
the momentum transfer and the hadron momentum. All the involved quantities have been 
conventionally defined in Ref. According to the latter convention, the Sivers function, 

f^j!(x,kx) |8j], and the Boer- Mulders function, h^ q (x,kT) [2], are formally defined through the 
following expression^]: 
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taking the proton polarized along the y axis; and 
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where S is the spin of the target hadron. The normalization of the covariant spin vector is 
S 2 = — 1, M is the target mass, ip q (£) is the quark field and Cg T is the gauge linko The gauge 
link contains a scaling contribution which makes the T-odd TMDs non vanishing in the Bjorken 
limit, as it has been shown in Refs. |12} \13\ I14j. 

The difference between the two functions is clearly seen comparing Eq. (pQ) and Eq. ([2]). 
The BM function counts transversely polarized quarks, hence the Dirac operator 7 + 7 2 7s in Eq. 
([2]), in an unpolarized proton. On the other hand, the Sivers function counts the unpolarized 
quarks, hence the Dirac operator 7 + in Eq. ([1]), in a transversly polarized proton, denoted by 
the transverse component S y in the proton state in Eq. (pQ). If there were no scaling contribution 
of the gauge link, the two T-odd functions, f^(x, kr) and h^ q (x, kr), would be identically zero. 

2. Quark model calculations 

Several quark model calculations of the Sivers and BM functions have been performed in the 
past years (see, i.e., Refs. [151 13 El US] ) and their phenomenology keeps attracting interest, 
as it is demonstrated by very recent important contributions [T9 l \20 \ [21] . 

The general formalism for the evaluation of time-reversal odd TMDs in quark models, 
presented in Ref. [221 123] 24], is reviewed here below. It is based on an impulse approximation 
analysis, where the Final State Interactions are introduced through a One Gluon Exchange 
mechanism, as depicted in Fig. [TJ The used quark models, i.e., the simplest version of the MIT 
bag model [25], and a non relativistic constituent quark model (NRCQM), do not contain explicit 
gluonic degrees of freedom. In the Sivers and BM functions the Dirac operators determine the 

1 a± = (ao ± o 3 )/a/2- 

2 The gauge link is defined as (00, £ - ) = Vexp (^—ig A + (r]~ , £t) drf J , where g is the strong coupling 
constant. This definition holds in covariant (non singular) gauges, and in SIDIS processes, as the definition of 
T-odd TMDs is process dependent. 




Figure 1. The process taken into account, together with its h.c, in the present evaluation of 
the T-odd TMDs. 



spin structure of the interaction, V mim2rn3rn4 between the quarks with initial spins ni2 and 
and final spins m\ and 7713, respectively. For instance, in the MIT bag model [25] . one gets, for 
the Sivers [BM] function 
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The imaginary part is taken only for the Sivers function calculation, performed changing the 
transversely polarized nucleon states from a transversity to a helicity basis. The coupling of the 
spins at the quark level appears clearly from the expressions 
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whose calculation has been performed assuming SU(6) symmetry for the proton state, for 
q = u,d. The interaction V(k,fa,qT) has been evaluated using the properly normalized fields 
for the quark in the bag [25], given in terms of the quark wave function in momentum space, 
which reads 



<Pm(k) 
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with the normalization factor N given in Ref. [25] The interaction is then 

V(k, fa, gT)m 1 m2,m 3 m4 = ~n <Pm x $ ~ <Zr) ^^ + T m ip m2 (k)^ {fa) 7°7 + <p m4 (fa - qr) , (6) 



where = 1 or 7 2 7s for, respectively, the Sivers and the Boer-Mulders function. The 

expressions in the NRCQM are similar to the one in the bag model. They read, being ^ the 



intrinsic proton wave function, 
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where ^ 5 = J2s z =-i 1 or =-1 1 &y f° r the ^M an d the Sivers function, respectively. The 
interaction here reads 

V(ki,k 3 ,q) mirn2m3mi = ^ w mi (A-i)7 + r /[hl u m2 (ki + q)u m3 (k 3 ) 7+ u m4 (fc 3 - g) , (8) 

with ii(/c) the four-spinor of the free quark states, arising from the Impulse Approximation 
analysis 

In Refs. [221 E3J 121], the T-odd TMDs have been evaluated in both the MIT bag model and 
the NRCQM assuming initially an SU(6) symmetry for the proton. 

In the former case, the qualitative results are understood calculating the coefficients (jU), in 
terms of which it is possible to realize what happens to the quark spins in the FSI process in a 
perfectly transparent way. In the latter case, one first has to re-express Eq. ([7]) in terms of the 
proton state. In spectroscopic notation and with the Jacobi coordinates, one has, using SU(6) 
symmetry 
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where \x)s 1S the standard SU(Q) vector describing the spin-flavor structure of the proton. In 
the bag model calculation, a non-zero Sivers function arises through the interference between the 
upper and lower components of the wave function. Something similar happens in the CQM. This 
time the interference is between the upper and lower components of the free spinors appearing 
as a consequence of the use of the Impulse Approximation (see Eqs. ([6l [8]) for the interaction 
and Refs. [221 |23[ [23] for details). Two different spin combinations contribute to the Sivers 
function. One of them comes from the spin-flipping of the quark interacting with the photon, 
i.e. the Y term in the MIT bag calculation 

firfabr) = Y^° 2 J '^^lC q + Y(q ± ,k T ) + C+-U(q ± ,k T )] 1 (10) 
with C a normalization factor, C q a weighting spin-flavor-color factor resulting from the matrix 
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elements @ and where Y/U(g±,hr) include the momentum dependent part] 

On the other hand, there are more spin combinations for the BM function|f| The first reason 
is that both non-flipping and double- flipping terms are important. The second reason is the sum 



3 See Eqs. (8-9) of Ref. [2 

4 See Eq. (13) of Ref. [24] 
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Figure 2. Comparison of the moments, f (h)f^ u ^ d \x) , Eq. (jlip . of the Sivers (red) and Boer- 
Mulders (blue) functions in the SU(6) NR Constituent Quark Model (left) and the MIT bag 
model (right). Short-dashed curves represent the d-quark distributions; full curves the w-quark 
ones. 



over the two spin states, i.e. S z = —1, 1. Due to the spin-flavor-color coefficients, i.e., due to the 
SU (6) symmetry assumption, the non-flipping term is more important than the double-flipping 
contribution. In effect, the latters are governed by the product of the two lower components of 
the bag wave function which encodes the most relativistic contribution arising in the MIT bag 
model. They turn out to be a few orders of magnitude smaller than the dominant ones, arising 
from the interference between the upper and lower parts of the bag wave function. This also 
happens if a proper non relativistic reduction of the gauge link, suitable for CQM calculations, 
is performed, justifying then the non relativistic approximation. 

In Fig. 2 the "first moments" of the Sivers and the Boer-Mulders functions, i.e. the quantities 

f# 1)Q to = J d 2 k T ^f^(x, k T ) , (11) 
are shown for u and d quarks in both the CQM and MIT bag model, using a a (^g)/(47r) — 0.13 

m- 

In Fig. 3 it is shown, in the MIT bag model case, the effect of neglecting the possible 
contribution of a spin-flip of the spectator quark (the one with initial spin projection in Fig. 
1), as was done in a previous calculation [15J. The proper inclusion of this term is crucial for 
the fulfillment of the Burkart Sum Rule, as it is explained here below. 

The above-described formalism has been easily extended, in the NR case, to the model of 
Isgur and Karl [26] ; allowing for a weak SU(6) breaking. The detailed procedure and the final 
expressions of the Sivers function in this model can be found in Ref. [22] . To evaluate numerically 
the Sivers and BM functions, the strong coupling constant g (i.e. a s (Q 2 )) has to be fixed. The 




h igure 3. The quantity f^ 1)uid) {x), Eq. (HU). Dashed curve: the results in the MIT bag 
model at //q, considering that a spin- flip can occurr only for the interacting quark. Full curve: 
the same, evaluated properly taking into account also the possible spin-flip of the spectator 
quark. 
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b igure 4. Left (right): the quantity /^(i), Eq. (JUJ). Dashed curve: the result in the 
Isgur and Karl model at /Xq. Full curve: the evolved distribution at NLO. Patterned area: the 
1 — a range of the best fit of the HERMES data proposed in Ref. |28| . 



prescription of Ref. |27j. is used to fix /j,q, according to the amount of momentum carried 
by the valence quarks in the model. Here, assuming that all the gluons and sea pairs in the 
proton are produced perturbatively according to NLO evolution equations, in order to have 
~ 55% of the momentum carried by the valence quarks at a scale of 0.34 GeV 2 one finds that 
Hi ~ 0.1 GeV 2 if Ag£g ~ 0.24 GeV. This yields a,(^)/(4vr) ~ 0.13 [27]. The results of the 
present approach for the first moments of the Sivers function are given by the dashed curves 
in Fig. 4, where they are compared with a parameterization of the HERMES data, taken at 
Q 2 = 2.5 GeV 2 . The patterned area represents the 1 — a range of the best fit proposed in Ref. 
|28j . The results are close, in magnitude, to the data, although they have a different shape: the 
maximum (minimum) is predicted at larger values of x. Actually /x 2 , is much lower, Q 2 = 2.5 
GeV 2 . For a proper comparison, the QCD evolution from the model scale to the experimental 
one would be necessary. Unfortunately, the Sivers function is a TMD PDs and the evolution 
of this class of functions is, to a large extent, unknown. In order to have an indication of the 
effect of the evolution, a NLO evolution of the model results has been performed, assuming, 
for f$ 1)Q (x), the same anomalous dimensions of the unpolarized PDFs. From the final result 
(full curve in Fig. 2), one can see that the agreement with data improves dramatically and the 
trend is reasonably reproduced at least for x > 0.2. Although the performed evolution is not 
exact, the procedure highlights the necessity of evolving the model results to the experiment 
scale and it suggests that the present results could be consistent with data, still affected by large 
errors. A discussion concerning the use of perturbative methods to relate the scale of model 
calculations to the scale of experimental data, and its impact on the calculation of T-odd TMDs, 
has been recently presented in Ref. |29| The deep understanding of the QCD evolution of TMDs 
is progressing fast [301 13T1 1321 1531 153] . 

Let us see now how the results of the calculation compare with the Burkardt sum rule [35] , 
which follows from general principles. The Burkardt Sum Rule (BSR) states that, for a proton 
polarized in the positive y direction, Y1q=u d(^) = with 

<A£) =~f Q dxJ dk T ^ft T Q (x, hr) , (12) 

and must be satisfied at any scale. Within our scheme, at the scale of the model, it is found 
(fcjf) = 10.85 MeV, (k£) = —11.25 MeV and, in order to have an estimate of the quality of the 
agreement of our results with the sum rule, we define the ratio r = \(k~) + — (k%}\ 

obtaining r ~ 0.02, so that one can say that the calculation fulfills the BSR to a precision 
of a few percent. One should notice that the agreement which is found is better than that 
found in previous model calculations, where the Sivers function for the flavor u was found to be 
proportional to that for the flavor d. For the MIT bag model, where the translational invariance 
is broken and the wave functions are not exact momentum eigenstates, the situatuion is slightly 
worse and the Sum Rule turns out to be violated of ~ 5 % (one should realize anyway that, in 
a previous calculation without spin flip of the spectator quark, the violation was found to be ~ 
60 %). 

The recent calculation of Ref. [18] , using overlaps of Light Cone Wave Functions in the Light 
Front of Dynamics, where the number of particles and the momentum sum rules are satisfied at 
the same time, fulfills the Burkart Sume Rule exactly, showing the importance of using exact 
momentum eigenstates in checking the transverse momentum of the quarks. 

3. The Sivers function from neutron ( 3 He) targets 

The experimental scenario which arises from the analysis of SIDIS off transversely polarized 
proton and deuteron targets is puzzling [36, 37]. With the aim at extracting the neutron 
information to shed some light on the problem, a measurement of SIDIS off transversely polarized 
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Figure 5. Left (right): the model neutron Collins (Sivers) asymmetry for 7r~ production (full) 
in JLab kinematics, and the one extracted from the full calculation taking into account the 
p p (dashed), or neglecting it (dotted). The results are shown for z=0A5 and Q 2 = 2.2 GeV 2 , 
typical values in the kinematics of the JLab experiments. 
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Figure 6. Left (right) The same as in Fig. 
the JLab experiments. 
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5, but at z=0.6, the maximum value accessed in 



He has been addressed [38J, and an experiment, planned to measure azimuthal asymmetries 
in the production of leading ir^ from transversely polarized 3 He, has been just completed at 
Jefferson Lab, with a beam energy of 6 GeV |39j. Another experiment will be performed after 
the 12 GeV upgrade of JLab |40| I41]. Here, a realistic analysis of SIDIS off transversely polarized 
3 He, presented in Ref. [12], is sommarized. The formal expressions of the Collins and Sivers 
contributions to the azimuthal Single Spin Asymmetry (SSA) for the production of leading pions 
off 3 He have been derived, in impulse approximation (IA), including also the initial transverse 
momentum of the struck quark. The final equations are rather involved and they are not 
reported here. They can be found in [32]. The same quantities have been then evaluated 
in the kinematics of the planned JLab experiments. Wave functions [43] obtained within the 
AV18 interaction [44] have been used for a realistic description of the nuclear dynamics, using 
overlap integrals evaluated in Ref. [IS], and the nucleon structure has been described by proper 
parameterizations of data [46J or suitable model calculations [47] , The crucial issue of extracting 
the neutron information from 3 He data will be now thoroughly discussed. As a matter of facts, 
a model independent procedure, based on the realistic evaluation of the proton and neutron 
polarizations in 3 He [48] . called respectively p p and p n in the following, is widely used in inclusive 
DIS to take into account effectively the momentum and energy distributions of the polarized 
bound nucleons in 3 He. It is found that the same extraction technique can be applied also in 
the kinematics of the proposed experiments, although fragmentation functions, not only parton 
distributions, are involved, as it can be seen in Figs. 5 and 6. In these figures, the free neutron 
asymmetry used as a model in the calculation, given by a full line, is compared with two other 
quantities. One of them is: 

where i stands for "Collins" or "Sivers" , A^ p ' is the result of the full calculation, simulating 
data, and d n is the neutron dilution factor. The latter quantity is defined as follows, for a 
neutron n (proton p) in 3 He: 

E ? e 2 q f q ' n(p) (x) D q ' h (z) 
dn{p){X ' Z) = Z N = P ,nZ q etfi> N (x)D^(z) (M) 

and, depending on the standard parton distributions, f q ' N (x), and fragmentation functions, 
D q,h (z), it is experimentally known (see [12J for details). A l n is given by the dotted curve in the 
figures. The third curve, the dashed one, is given by 

Ai ~ -L- - 2p p d p A^A , (15) 

Pn "n ^ ' 

i.e. 3 He is treated as a nucleus where the effects of its complicated spin structure, leading to a 
depolarization of the bound neutron, together with the ones of Fermi motion and binding, can 
be taken care of by parameterizing the nucleon effective polarizations, p p and p n . One should 
realize that Eq. (|13|) is the relation which should hold between the 3 He and the neutron SSAs 
if there were no nuclear effects, i.e. the 3 He nucleus were a system of free nucleons in a pure S 
wave. In fact, Eq. (|13p can be obtained from Eq. (|15p by imposing p n = 1 and p p = 0. It is 
clear from the figures that the difference between the full and dotted curves, showing the amount 
of nuclear effects, is sizable, being around 10 - 15 % for any experimentally relevant x and z, 
while the difference between the dashed and full curves reduces drastically to a few percent, 
showing that the extraction scheme Eq. (|15p takes safely into account the spin structure of 3 He, 
together with Fermi motion and binding effects. This important result is due to the peculiar 
kinematics of the JLab experiments, which helps in two ways. First of all, to favor pions from 



current fragmentation, z has been chosen in the range 0.45 < z < 0.6, which means that only 
high-energy pions are observed. Secondly, the pions are detected in a narrow cone around the 
direction of the momentum transfer. As it is explained in |42| . this makes nuclear effects in the 
fragmentation functions rather small. The leading nuclear effects are then the ones affecting 
the parton distributions, already found in inclusive DIS, and can be taken into account in the 
usual way, i.e., using Eq. ()15p for the extraction of the neutron information. In the figures, 
one should not take the shape and size of the asymmetries too seriously, being the obtained 
quantities strongly dependent on the models chosen for the unknown distributions [37] • O ne 
should instead consider the difference between the curves, a model independent feature which 
is the most relevant outcome of the present investigation. The main conclusion is that Eq. (|15p 
will be a valuable tool for the data analysis of the experiments [39} I40|. 

While the analysis of Ref. [32] has been performed assuming the experimental set-up of the 
experiment described in Ref. [39], but using DIS kinematics, a further analysis is being carried 
on [49] . to investigate possible nuclear effects related to the finite values of the momentum and 
energy transfers, Q 2 and u, in the actual experiment. In Ref. |49j . the description of Ref. [32] 
will be also improved, implementing a relativistic Light Front treatment to evaluate the nuclear 
polarized spectral function. Besides, the problem of possible effects beyond IA, such as final 
state interactions, will be addressed, and more realistic models of the nucleon structure, able to 
predict reasonable figures for the experiments, will be included in the general scheme. 
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